Multifractality and intermediate statistics in quantum maps 
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long-range statistics. Thus our results indicate with some 
generality the existence of a link between the statistics of 
eigenvalues and the multifractal properties of the eigen- 
functions. 

Let us start with the classical map defined on the 2- 
torus by <I> 7 : p = p + 7 (mod 1) ; q = q + 2p (mod 1), 
where (p, q), the coordinates in phase space, are the con- 
jugated action and angle variables and the bars denote 
the resulting variables after one iteration of the map. The 
quantization of this map yields a unitary evolution oper- 
ator acting on a Hilbert space of dimension N = \/(2itK) 
which can be expressed in momentum space by the N x N 
matrix @, % 
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We study multifractal properties of wave functions for a one-parameter family of quantum maps 
displaying the whole range of spectral statistics intermediate between integrable and chaotic statis- 
tics. We perform extensive numerical computations and provide analytical arguments showing that 
the generalized fractal dimensions are directly related to the parameter of the underlying classical 
map, and thus to other properties such as spectral statistics. Our results could be relevant for 
Anderson and quantum Hall transitions, where wave functions also show multifractality. 
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Statistical properties of wave functions and energy lev- 
els of quantum systems have been an important topic of 
research in the past decades with many applications to 
physical systems. It is now well-known that energy levels 
of e.g. quantum systems whose classical limit is chaotic, 
or disordered systems when eigenstates are extended, fol- 
low Random Matrix Theory (RMT). In this case, wave 
functions are typically ergodic and level statistics show 
level repulsion at short distances. Conversely, systems 
whose classical limit is integrable, or disordered systems 
in a regime of Anderson localization, show Poisson statis- 
tics of energy levels (without level repulsion), and wave 
functions are typically localized in phase space 

It has been realized only later that another universality 
class exists which is intermediate between the latter two. 
It can be observed in disordered systems at the Anderson 
transition or in certain systems whose classical limit 
is pseudo-integrable In this case, level statistics fol- 
low specific laws called semi-Poisson statistics, and wave 
functions generally show multifractal properties. This 
multifractal behavior has been extensively studied in the 
case of the Anderson transition 0, SB]) and has been 
also seen in quantum Hall transitions [7[. 

Recently, a simple model for intermediate systems was 
introduced which corresponds to a quantization of a cer- 
tain interval-exchange map Q. The model, although 
very simple, can display the whole range of semi-Poisson 
statistics when a parameter is changed. Moreover, a cer- 
tain randomization of this system was shown to yield a 
new model of Random Matrices with intermediate statis- 
tics @. 

Here, we examine multifractal properties of eigenfunc- 
tions for the Random Matrix model corresponding to 
intermediate quantum maps. We compute the inverse 
participation ratios (IPR), fractal dimensions, and sin- 
gularity spectra in a variety of regimes with different 
numerical methods. Using extensive numerical studies 
and analytical arguments, we show that the parameter 
of the model can be related to the fractal dimensions of 
the eigenfunctions, as well as to the spectral statistics. 
This Random Matrix model is known to span the whole 
range of semi-Poisson statistics for both short-range and 
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with <ftp = —2irp 2 /N. From this quantized map one can 
construct an ensemble of random matrices, taking <ft p as 
random variables uniformly distributed in [0, 27r[ [9(. The 
statistical properties of the pseudo-spectrum (the set of 
eigenphases) of U strongly depend on the value of the 
parameter 7. On the one hand, for generic irrational 
7, the spectral statistics of U are expected to follow 
those of the Circular Unitary Ensemble (CUE) of RMT 
if the (ftp are independent (non-symmetric case), or the 
Circular Orthogonal Ensemble (COE) if one imposes a 
symmetry <pN-p — 4>p- O n the other hand, for rational 
a/b, a variety of different behaviors are observed 
It was shown in Q that for aN = ±1 mod b the 
spectral statistics is of semi-Poisson type. In particu- 
lar the nearest-neighbor spacing distribution is given by 
P/s(s) = Aps@ e~^ +1 ^ s with parameter (3 = b — 1 in the 
non-symmetric case (/3 = 6/2— 1 in the symmetric case). 
For aN 7^ ±1 mod b, P(s) is still of intermediate type 
but with more complicated formulas [lol ] . Finally when 7 
is an integer the eigenphases are equally spaced and the 
spectrum is totally rigid. Thus the set of quantum maps 
U with rational 7 gives a random matrix ensemble with 
intermediate statistics (ISRM) whose spectral statistics 
correspond to natural intermediate distributions between 
Poisson and RMT, controlled by the value of 7. 

Multifractality properties of wave functions are de- 



2 



scribed by a whole set of generalized fractal dimensions 



D q . For a vector 



in an A^-dimensional 



Hilbert space, the multifractal exponents D q are defined 
through the scaling of the moments 
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The fractal dimension for q = corresponds to the di- 
mension of the support of the measure, here Dq — 1. The 
fractal exponent L>2 describes the large-size behavior of 
the IPR £ = l/J2iLi IV'il 4 ! which measures the extension 
of the state over the basis vectors. The multifractal 
exponents describe the behavior of the partition function 

N b 

Z(q,L) = J2ML) q <*L T «, r q = D q (q-l), (3) 
fc=i 

where the vector \ip) is divided into Nb = N/L boxes 
B k of size L, and fj, k (L) = J2ieB k IVM 2 , 1 < k < N b . 
The multifractal properties are alternatively character- 
ized by the singularity spectrum /(a), which is the frac- 
tal dimension of the set of points whose singularity ex- 
ponent is a. It is related to the function r q by a Leg- 
endre transform. Introducing the normalized measures 
(j,k(q,L) = Hk(L) q '/ 'J2i t 1 i( L ) q ') tne singularity exponent 
and the associated fractal dimension can respectively be 
obtained by fill ] 
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J2i Vi(q,L) log m(q,L) 



log(L/N) 



(4) 

Let us consider an ensemble of matrices of type ([1]) 
with rational 7 = a/6, in the non-symmetric case where 
all 4> p are independent. The mean IPR for eigenvectors of 
these matrices in p representation for different values of 
7 with denominator b = 7 is displayed in Fig.[T]as a func- 
tion of the matrix size. The IPR scales as N° 2 provided 
data corresponding to different values of aN mod b be 
treated separately. Indeed, when different matrix sizes 
are grouped into families, the results yield a linear be- 
havior of log(£) as a function of log TV, with the same 
slope £>2 for each family. More generally we observed 
that the fractal exponents D q are well defined if data are 
organized into families, and that they only depend on the 
denominator b of 7. 

We now proceed to compute the fractal exponents D q . 
A few of these exponents have already been computed in 
0, Qjj for the case aN = ±1 mod b. Here our aim is to 
characterize D q as a function of q. The quantities D q and 
/(a) are known to be difficult to compute numerically, 
especially for large q or a. In this work, we resorted to 
several different methods as a consistency check. We first 
opted for the usual method of moments. We computed 
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FIG. 1: (Color online) Mean IPR of eigenvectors of ([T]) as a 
function of the vector size N for 7 = 0/6 with 6 = 7 and 
a = 1 (filled symbols), a = 2 (half- filled squares), a = 3 
(empty squares). Straight lines correspond to the best linear 
fits. Logarithm is natural. 



average values of the moments ([2]) for different system 
sizes N ranging from ~ 2000 up to ~ 12000 to get rid 
as much as possible of finite size effects [13j . The fractal 
dimensions are extracted from the slopes of the graphs 
of log(^j IV^I 29 ) versus log N. Here the average is taken 
over all eigenvectors and random realizations of U (from 
200 realizations for N - 2000 to 1 for N ~ 12000). We 
also opted for the so-called canonical method [11] allow- 
ing to determine the f(a) spectrum directly from Eq. ([?]). 
For this method, the numerical computations were done 
on a single realization of size A^ ~ 13000, and 20 box 
sizes ranging from L = 10 to L = 0.1 AT (again different 
families of box sizes were treated separately). We also 
considered other approaches, such as the box counting 
method based on Eq. ([3]) ; they all give results intermedi- 
ate between the two previous methods. 

The results for D q are displayed in Fig. O For increas- 
ing b, the curve for D q tends slowly to the limiting curve 
where D q — 1 for all q, which corresponds to non-fractal 
wave functions. This is in agreement with the fact that 
7 tends for b — > 00 to an irrational number for which the 
system is expected to follow RMT. Figure [5] shows that 
D q is roughly linear in a relatively large interval around 
q = 0, and tends to limiting values D± oc for large \q\. 
The slope of D q at q = is displayed in the inset of 
Fig. [3] as a function of b. We found that the value of this 
slope is very accurately given by —1/b. Since Dq = 1, 
the first-order expansion of D q around q = reads 
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This expansion turns out to be valid in a quite large 
interval of q, whose size increases with b. As an example, 
the numerical values of D\ and D2 together with the 
linear approximation ([5]) are shown in Fig. [3J When b 




FIG. 2: (Color online) Fractal dimension D q computed with 
the canonical method for 7 = 1/6 with b — 2 (blue circles), 
6 = 3 (red squares), b — 10 (green triangles). Solid lines show 
the linear approximation D q = 1 — q/b. Inset: anomalous 
exponent A q as a function of q (solid line), together with its 
symmetry with respect to q = 1/2 (dashed line) for 6 = 2. 

is increased, Eq. ([5]) is verified with higher and higher 
accuracy. Of course Eq. breaks down for large \q\ 
since D q is bounded. 

To get an understanding of why formula ^ holds, 
we note that general arguments for critical systems pre- 
dict that the multifractal properties of the eigenstates for 
q = 2 are linked to the spectral statistics through a rela- 
tion between the correlation dimension D2 and the level 
compressibility x 0]> 



Numerical results for the power-law random banded ma- 
trix ensemble [5| have revealed that this relation is ex- 
tremely well verified in the regime of weak multifractal- 
ity (large bands): in this model, the fractal dimension 
evolves linearly with respect to q as D q — 1 — nq, where 
K is inversely proportional to the width of the central 
band (and in this particular case can be also related to 
the level compressibility). However, for smaller bands 
Eq. ([6]) was clearly violated. Suppose Eq. ((6]) holds in our 
case. For ISRM (Eq. (JlJ), the level compressibility can 
be estimated analytically. It is given by the value of the 
two-point correlation form factor ^(t) = | Tr U n \ 2 /N, 
with t = n/N, for n/N -> 0. Following ||, we note that 
in the semiclassical limit N — > 00 and fixed n, the trace 
TrU n is asymptotically equal to TrV n , where V n is the 
quantization of the nth iterate of the classical map, and 

^ N-l N-l 

Tr V n = — exp(m0 p ) ^ exp(2i7m7fc). (7) 

p=0 fc=0 

The modulus squared of the first sum yields ss N when 
all <j) p are random. The second sum is a geometric sum: 



FIG. 3: (Color online) Information dimension Di (blue cir- 
cles) and correlation dimension D2 (red diamonds) as a func- 
tion of the denominator of 7 = 1 /b. Full (empty) symbols are 
numerical values obtained by the method of moments (the 
canonical method). Solid and dashed lines are the theoretical 
curves 1 — 1/6 and 1 — 2/6 for Di and D2 respectively. The 
values obtained by e.g. the box counting method lie in the 
shaded domain in between. Inset: slope of D q at the origin 
q = (triangles) and curve —1/6 (solid line). 

for 7 = a I b, it is equal to N if n is divisible by b and 
to O(l) otherwise. Thus, K2{n/N) ~ 1 provided n is 
divisible by b, K2(n/N) ~ in all other cases. The level 
compressibility is then given by the time averaged form 
factor 

1 " 1 

x = K 2 {0)=\im lim - V K 2 (n'/N) r» -. (8) 

n^oo N^oo n * — ' 
n' = l 

Inserting this value of \ m to Eq. (J6j) we get D2 ~ 1 — 2/6, 
which corresponds to Eq. (|5|) for q — 2. A simple linear 
interpolation between this value for D2 and Dq = 1 yields 
Eq. $5§. We note that for small b (strong multifractality) 
Eq. ([6]) breaks down but Eq. ([5]) is still valid for smaller 
q values. 

Before moving to the study of the singularity spec- 
trum, we briefly discuss symmetry properties of D q . It 
was suggested in [6( that the anomalous exponents A q , 
defined by A q = (D q — l)(q — 1), approximately follow 
the symmetry relation A q = Ai— q . This was shown to 
hold for the Anderson model with good accuracy over a 
large interval of q values. It is not the case in our system. 
As an example, the inset in Fig. Ogives A g and Ai_ 9 for 
7 = 1/2. For values of q where the exponents D q have a 
linear behavior the symmetry relation holds, as it should 
since any linear D q necessarily fulfills it. Outside the 
linear regime, the relation is not verified anymore. 

We now turn to the singularity spectrum /(a). For 
7 = a/b, the expression obtained using Eq. ([5]) is 

/(a)«l-j(a-l-i) • (9) 
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FIG. 4: (Color online) Singularity spectrum /(a) for 7 = 1/6 
with 6 = 2 (blue circles), 6 = 3 (red squares), 6 = 10 
(green triangles). Solid lines show the parabola Q. In- 
set: singularity spectrum width w = Q max — a m i„ as a 
function of 6. The dashed line shows the best linear fit 
logw = —0.530 log 6 + 1.247. Logarithm is natural. 

It reaches its maximum at a(q = 0) = 1 + 1/6. Since 
Eq. is valid around q — 0, we expect Eq. §§§ to be 
accurate around a(0). Figure 0] shows the singularity 
spectrum, numerically computed using Eq. Q, together 
with the theoretical estimate Eq. ([9]). The data dis- 
played show that ([9|) approximates the singularity spec- 
trum with good accuracy over a large interval of values 
of a around a(0). As 6 — > 00, the curve for /(a) gets 
closer and closer to a single point at a = 1, once again 
corresponding to the non- fractal case of RMT. 

We previously showed that the theory ([5]) accurately 
describes the behavior of the moments in Eq. j2|) for small 
q. For large values of g, this is bound to break down since 
D q converges to a finite asymptotic value for q — > ±00. 
Similarly /(a) should have vertical asymptotes at some 
limiting values a max and a m i n , while ([9]) is the equa- 
tion of a parabola. However numerical data are not far 
from the theory ((9J, and some of the features of /(a) are 
well captured by this estimate. For example, the inset 
of Fig. 2] shows that the width w = a max — a m i n of the 
singularity spectrum scales as ~ l/6 a53 (best fit). This 
is close to the scaling 1 / y/b of the difference between the 
two intersections of the parabola (J9j) with a straight line. 

We finish by noting that in the symmetric case 
= <Pp m Eq. |l|) we have performed similar 
computations, getting very similar data. In particular, 
Eqs. ([5]) and ([9]) are valid in this case as well. 

In conclusion, we have studied multifractal properties 
of eigenfunctions for intermediate quantum maps. Al- 
though data corresponding to system sizes N with dif- 
ferent values of aN mod b should be treated separately, 



they give the same value for D q and f{a). Our results 
show that for an interval of q values whose size increases 
with b, the fractal exponents can be related explicitly to 
the parameter 7 = a/6 of the map through Eq. ([5]), and 
thus to spectral statistics. A similar result holds for the 
singularity spectrum through Eq. Thus in such a 

system, fractal exponents and singularity spectrum are 
related to the spectral properties over a wide range of 
fractal dimensions. Interestingly enough, our relation is 
still valid for small q even when Eq. (J6j> for q = 2 does not 
hold. As our system corresponds to a Random Matrix 
model covering the whole range of semi-Poisson statis- 
tics, both at short-range and long-range, we can expect 
our results to display some generality. It will be interest- 
ing to study if similar results apply to other intermediate 
systems, and other physical systems where wave func- 
tions are multifractal, such as condensed-matter systems 
at the Anderson or quantum Hall transitions. 
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